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Abstract: We compute the one- loop effective potential for the Polyakov loop on S 3 x S 1 
for an asymptotically free gauge theory of arbitrary group G and a generic matter content. 
We apply this result to study the phase structures of G2, SO(N) and Sp(2N) gauge theories 
which turn out to be in qualitative agreement with the results of lattice calculations. On 
x S , at zero temperature, the Polyakov loop is zero for kinematical reasons. For 
small but non-zero temperature, the Polyakov loop is still zero if the gauge theory has 
an unbroken center, while it acquires a small vacuum expectation value for gauge theories 
whose center is trivial or explicitly broken by the presence of dynamical matter fields. At 
high temperatures, the saddle point structure of the effective potential is different from 
the low temperature case suggesting that the theory is in the deconfined phase. At finite 
N, the Polyakov loop is non-zero in the high temperature phase only in theories with no 
unbroken center symmetry, consistent with screening of the external charge introduced by 
the Polyakov loop and Gauss law on a compact space. 
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1. Introduction 

Quantum Chromodynamics (QCD) is a confining theory: the elementary particles of the 
theory, quarks and gluons, never appear as final states of strong interactions. This property 
holds for temperatures below a critical value T c , above which the system is deconfined, and 
quarks and gluons form a plasma. Confinement is a property of both QCD and the pure 
gauge theory without any matter. To gain a better understanding of this phenomenon, it 
is useful to study confinement in pure Yang-Mills theory. For pure SU(N) gauge theory, 
the existence of a deconfinement phase transition has been rigorously proved in [1]. 

The finite temperature theory is formulated on a four dimensional Euclidean space 
with one compact dimension. The length f3 of the compact dimension is the inverse of the 
temperature T. In the confined phase, the free energy of an isolated static quark is infinite. 
The free energy F is related to the vacuum expectation value of the Polyakov loop 

L R (x) = ±-tr R (>/o /T ^o(W) , (1.1) 

where Aq is the vector potential in the compact direction (parameterized by the coordinate 
t), g is the coupling of the theory, d,R is the dimension of the irreducible representation 
under which the quark transforms 1 and x is the position of the quark. If we define the 
vacuum expectation value (vev) of L (we drop the representation index from now on) as 

< L > = ^(/ md 3 x^ , (1.2) 

1 In QCD, this is the fundamental representation; it is useful however to study a more general case, in 
which the quark is in a generic representation of the gauge group. 
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the relationship between L{x) and the free energy of the quark is given by [2] 

F = -T\og(L). (1.3) 

In SU(N) gauge theories, the behavior of (L) across the deconfinement phase transition 
can be explained in terms of a symmetry. This is the symmetry of the system under the 
Zjy center of the group, which leaves the Lagrangian unchanged. Under a transformation 
z £ 7*n, (L) — ► z(L). In the deconfined phase, F is finite, i.e. (L) ^ 0. This implies 
that in the deconfined phase the center symmetry is broken. Conversely, in the confined 
phase F is infinite, (L) = and the center symmetry is not broken. This leads to the 
natural interpretation that confinement reflects a change in the property of the vacuum 
under the Z at center symmetry, the vev of the Polyakov loop being the order parameter for 
the transition. This observation lead Svetitsky and Yaffe to the well-known conjecture [3] 
that relates the deconfinement phase transition in a four dimensional SU (N) gauge theory 
to the order-disorder transition in three dimensional Zjv spin systems: since the properties 
of the deconfinement transition are determined by the underlying Z n symmetry, the gauge 
and the corresponding spin system are in the same universality class when both transitions 
are second order. 

Since confinement is a non-perturbative phenomenon, a fully non-perturbative ap- 
proach such as putting the theory on a space-time lattice is mandatory for quantitative 
studies. For SU(N) gauge theory, recently it has been shown in [4] that important in- 
sights about the physics of color confinement and the phase structure of the theory can 
be obtained by studying the system on a 5 3 x S 1 manifold, where the radius (3 of the S 1 
is connected to the temperature in the usual way (0 = l/T) and the radius R of the S 3 
provides the IR cutoff at which the running of the gauge coupling freezes. When 1 / R is 
much larger than the dynamical scale of the theory, the coupling is small and perturba- 
tive calculations are reliable. It is then possible to compute an effective potential for the 
Polyakov loop at one loop order in perturbation theory. Evaluating the partition function 
of the system for j3 <C R (high temperature regime) and (3 S> R (low temperature regime) 
yields information about the phase structure of the theory. In confining gauge theories, 
for a finite number of degrees of freedom a cross-over between the confined and deconfined 
phase is observed, with the two phases being characterized by a different structure of the 
minima of the free energy (or, using a more field theoretical terminology, the effective po- 
tential). In order to verify the existence of a real phase transition at some critical value of 
the temperature, an infinite number of degrees of freedom is needed. In standard thermo- 
dynamical calculations, this is achieved by taking the infinite volume limit of the system. 
This limit is obviously not possible for S 3 x S 1 , since the calculations rely on R being small. 
Instead, the thermodynamic limit is achieved here in another way: by sending to infinity 
the number of elementary degrees of freedom of the theory (one example in SU(N) gauge 
theory is to take the limit of infinite number of colors N with fixed g 2 N [5]). Thus, the 
existence of a phase transition can be rigorously proved only in some large N limit; for 
finite N or for exceptional groups, weak coupling calculations on a 5 3 x S 1 manifold can 
only discover the existence of two regimes (confined and deconfined), but cannot establish 
whether they are separated by a cross-over or a real phase transition in infinite volume. 
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Weak coupling calculations on S 3 x S* 1 have been used for determining the phase 
structure of various super symmetric and non-supersymmetric gauge theories [6-10], treat- 
ing each case independently. In this paper, we shall show that a general expression can 
be derived from the representation theory of the Lie Algebra of the gauge group that al- 
lows one to compute the one-loop effective potential for a generic gauge theory containing 
bosons and fermions in arbitrary representations. Using our result, it is straightforward to 
work out the phase structure of non-Abelian gauge theories for any gauge group and given 
matter content. We shall show that our computation reproduces the previously known 
results and we will use it to determine the phase structure of other gauge theories that 
can help us to understand the physics of color confinement, like G2 gauge theory, with no 
center, and SO(N) gauge theory, whose center is either trivial (N odd) or Z2 (N even) 2 . 
While a priori it is not obvious that the small volume calculations reproduce features of 
the thermodynamical limit, comparisons with other techniques (such as the lattice, as we 
will show in this paper) provide evidence that this is indeed the case. 

Our calculation show that in general the minimum structure of the effective potential is 
different at low and high temperatures. This means that a cross-over (in the finite volume 
case) or possibly a phase transition (which can arise only if the thermodynamical limit 
can be taken, e.g. by sending to infinity the number of elementary degrees of freedom of 
the theory) separates those two regimes, which are naturally identified respectively with 
the confined and deconfined phase. However, in both regimes the details of the minimum 
structure depend on the center symmetry. A symmetry cannot be spontaneously broken 
on finite volume, so if there is a non-trivial center that is not broken explicitly by the 
matter content of the theory, the Polyakov loop cannot develop an expectation value at 
any temperature. Why the Polyakov loop vanishes in these theories differs in the low 
temperature and the high temperature phase. In the low temperature phase, the effective 
potential for the Polyakov loop will have only one saddle point in which the Polyakov 
loop will be zero. At high temperatures, there are multiple saddle points of the effective 
potential, each related by the action of the center. In each of the saddle points, the 
Polyakov loop is non-zero. However, since we are on finite volume, these saddle points are 
not superselection sectors and therefore we have to sum over all of them, which yields a zero 
value for the Polyakov loop. It is plausible to assume that the existence of the saddle points 
at high temperatures survives the infinite volume limit in which, we can restrict ourselves 
to one of the multiple saddle points, thereby breaking the center symmetry and signaling 
a confinement /deconfinment transition. The situation is different when the theory has no 
unbroken center subgroup. In that case, at both low and high temperatures, there is one 
global minimum of the effective potential. In the low temperature phase, the Polyakov loop 
is small {(D(e~PI R )) and is zero at zero temperature whereas in the high temperature phase, 
we find that it develops an 0(1) non-zero expectation value. This can be understood in 
terms of screening the external charge introduced by the Polyakov loop. For theories with 
a non-trivial center, the external charge cannot be screened by gluons and quarks popping 
out of the vacuum. However, if there is no center symmetry, any external charge can be 

2 In this paper, we will exclusively study SO(N) groups and not their covering groups, Spin(N). 
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screened by an appropriate number of gluons and quarks. 

The paper is organized as follows. In Sect. |2| we compute the effective potential for a 
generic Lie group on a S 3 x S 1 manifold and derive the behavior of the Polyakov loop at 
high and low temperature in terms of the root and weight vectors of the group and the 
matter content of the theory, emphasizing how our result reproduces previous calculations 
and how it can be used to analyze other gauge theories. In Sect. [3|, we specialize to the G2 
case, showing that the Polyakov loop is zero at sufficiently low temperature and different 
from zero at high temperatures. The same result is worked out for SO(N) gauge theories 
in Sect. ||. In Sect. || we show that our results are compatible with the non-Abelian Gauss 
law. Finally, Sect. || sums up our findings and their implications for our understanding of 
color confinement. 

2. Confinement on S 3 for an arbitrary gauge theory 

In the canonical ensemble, the thermal partition function of a quantum field theory is 
equal to the Euclidean path integral of the theory with a periodic time direction of period 
(5 = A and anti-periodic boundary conditions for the fermions along the temporal circle. 
We will consider a gauge theory with a semi-simple group G with matter in arbitrary 
representations on a spatial <S 3 . In the following, we will assume that the theories we want 
to study are asymptotically free. To understand the phase structure, we need to calculate 
the free energy which is just the log of the partition function. Therefore, we perform a 
Euclidean path integral for the theory on x S^. 

We take 1/R to be large compared to the dynamical scale of the theory , so the theory 
is at weak coupling. The S 1 topology allows non-trivial flat connections F^ u = that can 
be parameterized in a gauge invariant way by the vev of the Polyakov loop. The set of 
flat connections are the lightest degrees of freedom of the theory, and other excitations are 
separated in energy by a gap proportional to 1/R. In fact, the flat connection is a zero 
mode of the theory, in the sense that the quadratic part of the action does not depend on it. 
The low energy phase structure can be studied using a Wilsonian approach where heavier 
degrees of freedom are integrated out to give an effective action for the flat connection. In 
the weak coupling regime, we can perform a perturbative expansion in the coupling and the 
one loop result provides a good approximation to the effective potential. Our discussion 
will generalize previous results [7,11,12]. 

We can use time-dependent gauge transformations to fix Aq to be constant and in the 
Cartan subalgebra, 

r 

A = Y,H a C a , [H a ,H b } = 0\fa,b (2.1) 
a=l 

where r is the rank of the group. The vev of the Polyakov loop is 

L= tr exp (J A^j . (2.2) 

Notice that Aq is anti-Hermitian in this expression. A coupling constant factor has been 
absorbed in the field as well. This means that in the classical Lagrangian the fields will 
not be canonically normalized, but there is a 1/g 2 overall factor. 
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We will now consider adding matter fields to the gauge theory in arbitrary represen- 
tations that can be obtained as tensor products of the fundamental (and anti-fundamental 
in case it exists). The fundamental representation corresponds to a mapping of the group 
to the (complex) general linear group acting on a vector space of dimension F. So any 
representation can be mapped to a tensor (we are ignoring the spinor representations of 
SO(N) groups). The tensor can have covariant or contravariant indices, and the adjoint 
action can be defined as (summation understood) 

Ad(A) X = AfX^M +■■■ + A* n X h "i m - Al'x'-j" M 1 Xj 1 '"'™ . (2.3) 

For SU (N) groups, the Cartan subalgebra admits a diagonal form in terms of the funda- 
mental weights Pi, i = 1,...,F, where F is the dimension of the fundamental representa- 
tion. 

(2.4) 

The fundamental weights need not be linearly independent; this will happen only if F = r, 
the rank of the group. Real groups like SO(N) and Sp(N) cannot be diagonalized using 
gauge transformations. However, once we have fixed the flat connection to be constant, 
we can do a global similarity transformation to diagonalize it. Gauge invariant operators 
remain unchanged, so this transformation is allowed. In the case of SO(N) this is done by 
a U(N) transformation, while in the case of Sp(N), it is a U(2N) transformation. In the 
diagonal basis, the adjoint action simplifies to 

(n m \ 

a=l 6=1 / 

which is 

(n m \ 

XX " XX X ^ = [( '' • ^ X t't ■ (2-6) 
a=l 6=1 / 

Here, C , V and /x m ,n are r dimensional vectors and the dot product is the standard 
Euclidean product. The weights of any representation are linear combinations of the fun- 
damental weights, so fl m ,n = (X)o=i ^« — YlbLi i s a weight of the tensor representation 
with m covariant and n contra- variant indices. 

The covariant derivative acting on a field in the background field of the Polyakov loop 

is 

D^X = d^X - 8^kd(A) X. (2.7) 

To calculate the one-loop effective potential, we expand the gauge field, matter and the 
ghosts in harmonics on S 3 x S 1 and keeping only the quadratic terms, we integrate out the 
various modes. The one loop potential will be in general 

V ^ = ~\ E(" 1 ) FflTr s ;5 xSi log (-D 2 0R - A) , (2.8) 

R 

where R denotes the gauge group representation of the different fields we have integrated 
out, Fr is the fermionic number and the trace is over the physical degrees of freedom (zero 
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modes excluded). This can be written as 

= EW«-<?)' (2-9) 

R Pr 

where (1r are the weights of the representation and the potential Vr can depend on the 
Lorentz representation of the field. In flat space (T n x R d case) Vr is just an universal 
function times the number of physical polarizations (see [7], for instance). On S 3 the 
situation is different because the Kaluza-Klein decomposition depends on the spin of the 
field. The one-loop potential for several SU(N) theories on S 3 x S 1 has been computed 
previously [4,8]. We will generalize the computation to other groups and cast the results 
in a simple form. 
In general, 

Tr s3xSl log {-D 2 0R - A) = £ EE^) lo § U + (jT + /** • C ) 2 J ' ( 2 - 10 ) 

where Mr is the number of physical polarizations of the field and Et is the energy of the 
Kaluza-Klein mode ipi satisfying Aipi = —efipi (A is the Laplacian on 5 3 ). Mr and eg 
depend on the field type that is being integrated out in a way that we list below. 

(i) Scalars. Scalar field can be conformally or minimally coupled. The set of eigen- 
vectors is the same in both cases. However, their eigenvalues are different. For confor- 
mally coupled scalars 3 we have e% = R~ l {£ + 1) whereas for minimally coupled scalars 
El = R~ 1 y / £(£ + 2). In both cases, there is a degeneracy Mr{£) = (£ + l) 2 with £ > 0. 

(ii) Spinors. For 2-component complex spinors, we have En = R~ x (l + 1/2) and 
M R (£) = 2£(£ + 1/2) with £ > 0. 

(iii) Vectors. A vector field Vi can be decomposed into the image and the kernel of 
the covariant derivative: Vi = ViX + Bi, with V l Bi = 0. The eigenvectors for Bi, have 
e t = R~ 1 (£ + 1) and M R {£) = 2£{£ + 2) with £ > 0. On the other hand, the part V»x has 
E t = R~ 1 v / £(£ + 2) with degeneracy M R {£) = {£ + l) 2 but with £ > only. 

Notice that both spinors and vectors have no zero (£ = 0) modes on S 3 . This is why 
in a pure gauge theory the only field with a zero mode is Aq which is a scalar on S 3 . 

It is a standard calculation using Poisson resummation to show that up to an infinite 
additive constant, 

oo oo 1 

TWilog(-I>g fl - A) =EE A/ ^){^-E- e " n/te ' cos(nP fiR-C)} . (2.11) 

P.r e=o n=l 

The first term here involves the Casimir energy and since it is independent of C will play 
no role in our analysis and we will subsequently drop it. Ignoring this term, the potential 
Vr is given by 

00 sin 2 ( 



n 

n=l 



3 These have a mass term involving the Ricci scalar of the manifold, in this case simply R 1 . 
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The sum on n can be done explicitly and the result is, up to a constant 

v ^-^ me),o i i+ ^f)) <2 - i3) 

This is the contribution to the effective potential obtained by integrating out a particular 
field. The information about the type of field being integrated out (its Lorentz transfor- 
mation properties etc) is encoded in Mr{£) and £(_ as explained above. 

In the remainder of this section, we will restrict ourselves to a pure gauge theory with 
no matter fields. Then, the basic degree of freedom is the non-Abelian gauge potential 
which decomposes into a time component Aq which transforms as a minimally coupled 
(massless) scalar on S 3 and a vector Ai on S 3 . We can write A{ = VjC + Bi such that 
Vj-B* = 0. Gauge fixing (ViA 1 + DqAq^J = by the usual Fadeev-Popov procedure leads 
to the existence of ghosts. 

Then, integrating out the Aq, Bi, C and ghost modes can be performed using the 
following table: 

field S 3 KK-mode e t Mr{£) 
B % £>0 (£+l)/R 2£(£ + 2)x(+l) 

C £>0 yj£{£ + 2)/R (^ + l) 2 x(+l) 

ghosts £ > v / £(£ + 2)/R (£ + l) 2 x (-2) 
A £>0 y/£(£ + 2)/R (^ + l) 2 x(+l) 

Notice that for £ > 0, the modes C, ghost and Aq cancel and the only contributions is 
from divergenceless modes, Bi. The constant £ = modes for C, ghosts and Aq give the 
following contribution: 

Tr , . cos(nx) , 9 fx\ 

V{x)~-J2— log sin 2 (-J (2.14) 

n=l 

This contribution is precisely the logarithm of the Jacobian that converts the integrals over 
C into the invariant integral over the group element g = exp(j,p\^.C^j , yielding the correct 
gauge invariant measure in the path integral dfi(g), g £ G. 

dCl[ S in(^^-)...c fd^g)-... (2.15) 

MAd V / J 

The contribution from the divergenceless part of Ai which we have been denoting by Bi 
can be easily written down using ( 2.13| ) and reading off the corresponding and Nr{£) 



from the table above. Then, the total effective potential is 

^gaugc(C) = V (PAd • C) , (2.16) 



where 



MAd 



V(x) = -log S in(f)| + E/(/ + 2)log^l + s J 2 n2 ^ | . (2.17) 
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The two terms in the effective potential come with opposite signs: the first term corresponds 
to a log repulsion coming from the I = modes whereas the second term is an attractive 
term. It is the interplay between these two types of contributions that leads to an interesting 
phase structure for the theory. 

In the low-temperature limit (3/R^> 1, the attractive contribution to the potential is 
exponentially supressed 



V(x) ~ — log 



M „ , *) | + J2 41(1 + 2) sin 2 (|) e"^ + . . . . (2.18) 

l>0 



In the high-temperature limit fi/R <C 1, the lower mode terms become stronger and cancel 
the repulsive term except very close to the origin x ~ (3/R 



V(x) ~ - log sin - + 6 log sin - + . . . . (2.19) 



Therefore, we expect a change of behavior with increasing temperature. At low temper- 
atures the expectation value of the Polyakov loop is zero up to exponentially supressed 
contributions. This can be seen from the fact that the repulsive term becomes the group 
measure ( 2. 15) ) and, since the Polyakov loop corresponds to the character of the group 



associated to the representation R, 

(L) = \ J d»(g) XR (g) = . (2.20) 

In the zero temperature limit this is a kinematical statement: there cannot be unscreened 
charges in a compact space. Therefore, confinement is just a kinematical fact. The deriva- 
tion of the Haar measure in S 3 x S 1 was known for unitary groups [4], but our discussion 
is completely general and includes any semi-simple Lie group. 

It is interesting to notice that at higher temperatures the expectation value of the 
Polyakov loop could be non-zero even in finite volume. An intuitive explanation follows 
from considering the Polyakov loop as the insertion of an external source in the system. 
At zero temperature pair production from the vacuum is completely frozen out, so it will 
cost an 'infinite amount of energy' to introduce a charge in finite volume. However, at 
higher temperatures the system is in an excited state where the charges present in the 
medium can screen the external source. Whether this actually happens or not depends on 
the details of the group and the type of source we introduce. In gauge theories with unitary 
groups this issue boils down to whether the center of the group is broken or not. In pure 
Yang-Mills the center is unbroken, so the expectation value of the Polyakov loop should 
be zero at any temperature. On the other hand, in a theory with quarks the center is 
completely broken, so the Polyakov loop can be nonzero at high temperatures. In theories 
with an unbroken center, the Polyakov loop has to be zero due to kinematics. However, at 
high temperatures the main contributions to the path integral will come from a sum over 
saddle points corresponding to the minima of the potential. For very large temperatures 
the measure can be approximated by a sum of delta functions, so the value of the Polyakov 
loop is just the sum at the different minima. At each of these points the Polyakov loop 



-8- 



has a nonzero value and it is plausible to assume that there is a map between them and 
the possible vacua of the gauge theory in the deconfined phase in infinite volume, where 
kinematical constraints do not apply. In this sense, by restricting ourselves to only one 
saddle point at high temperatures thereby getting a non-zero value of the Polyakov loop, 
we can interpret these results as a confinement-deconfinement phase transition at weak 
coupling. 

Strictly speaking there is no transition at finite volume, in the high temperature phase 
we should sum over all the vacua so the average of the Polyakov loop will be still zero 
or we would go smoothly from a zero to a non-zero value as we increase the temperature. 
However, it is possible to show there is a sharp transition in the large ./V limit [4]. Of course 
this is not applicable to the exceptional groups but only to the classical groups. 

2.1 Large N phase transition 

In the ensuing analysis, it is convenient to write the potential in terms of the Polyakov loop 
eigenvalues Aj = v"i ■ C, where V{ are the fundamental weights. This immediately follows 
from the fact that the weights of any tensor representation are linear combinations of 
the fundamental weights, and the contributions from integrating out such representations 
to the effective potential will depend on corresponding linear combinations of \. In the 
N — > oo limit for the classical groups, we can approximate the discrete set of eigenvalues 
by a continuous variable 8 and substitute the sums over eigenvalues by integrals in 9 



The domain of integration is determined by the periodicity of the Polyakov loop eigenvalues. 
The potential becomes a functional of the eigenvalue distribution, and the system is well 
described by the saddle point that minimizes the potential. When the parameters of the 
theory change, such as the ratio /3/R, the dominating saddle point can change, leading to 
a phase transition. 

For pure U(N) Yang-Mills, the potential is given by equation fl2.17|) with x = /2 a dj-C = 
Aj — Xj. After converting the summation to the integral via the density of eigenvalues, we 
obtain 




(2.21) 



with an eigenvalue density p(9) > 0, such that 






where q = e@/ R and zy is the partition function of divergenceless vector fields in S 3 




(2.24) 
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We can expand the eigenvalue distribution in complex Fourier components p n 

1 ' ^ 

P(8) = 7T + Yj Re ^«) cos ( n9 ) + lm (Pn) sm{n6) (2.25) 
^ 1 1 

n=L n=l 

and perform the integrals. We are left with a quadratic potential for each Fourier compo- 
nent 

V = iVV f; {1 ~ Zviqn)) \p n \ 2 . (2.26) 

n=l 

For low enough temperatures, zy{q n ) < 1 for all n and all the coefficients in the summation 
are positive so the configuration of minimal energy is the uniform distribution p{0) = 
l/(27r). In this case the expectation value of the Polyakov loop vanishes, giving a weak 
coupling version of the confining phase. When we increase the temperature, the coefficient 
of the n = 1 component first becomes zero at a critical value q c = 2 + y/S and then negative 
beyond this point. If we continue increasing the temperature, eventually larger Fourier 

(n) 1/n 

components become unstable beyond the points qi = q c of marginal stability. At the 
critical value q c the eigenvalue distribution jumps to a qualitatively different saddle point 
characterized by the appearance of a gap in the eigenvalue distribution, i.e. there is a 
range of values on 6 over which that function p{9) has no support. At temperature just 
above q c , 

p(0) = ^(l + cos(0)) . (2.27) 

The expectation value of the Polyakov in this saddle point is no longer zero, so this tran- 
sition is a weak coupling analog of the confinement /deconfinement phase transition. 

2.2 Some examples 

We would like to emphasize that the formulas derived in this section are completely general 
on 5 3 x S 1 for any semi-simple gauge group and for fields in any representation. We can 
cast the complete one-loop potential as a sum of different terms determined by the spin 
of the field and the weight system of the group representation. Possible examples are 
pure gauge theories, with unitary, orthogonal symplectic or exceptional groups. We will 
comment some of them in the next sections. 

Here we give a small list of examples in terms of the eigenvalues of the Polyakov loop 
Aj = Vi ■ C for unitary groups. The phase diagram of most of them has been analyzed 
previously in the large- iV limit [6,9], showing an interesting structure in terms of the 
expectation values of the Polyakov loop. All of them exhibit a deconfinement transition at 
some critical temperature T c . We will use the potential generated by Weyl fermions 

The 7r/2 addition in the argument is for anti-periodic boundary conditions on fermions. 
Then, we can write the following potentials: 
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• QCD: SU(N) Yang-Mills with Nf massless flavors (or theories with fundamental 
matter in general) 

N 

V(X) =J2 v(\ i -X j )- 2Nf (Aj) • 

At high temperatures, the flavor contribution to the potential favors the Aj = 
minimum over other center-breaking minima. At low temperatures the flavor po- 
tential is exponentially supressed, and only the repulsive term from the pure gauge 
contribution is relevant. 

• Supersymmetric theories: M = 1 SU(N) (or theories with adjoint matter in general) 

Notice that when I — > oo, fermionic and bosonic contributions cancel at leading order 
if we choose periodic boundary conditions for fermions. At high temperatures all the 
center-breaking minima become deeper, so center symmetry is unbroken. At large 
N, the system will sit at one of these center-breaking minima. At low temperatures 
the fermionic contribution is exponentially supressed. 

• Orientifold theories [13]: SU(N) gauge theories with a Dirac fermion in the anti- 
symmetric (or symmetric) representation (or theories with two-index representations 
in general) 

v(x) = J2v(x i -x j )-J2v F (\ i + x j ) -U^VFixA . 

iy^j iy^j \ 1 = 1 / 

The term inside the parentheses only appears in the symmetric representation; in 
the large- N limit its contribution is sub-leading. At high temperatures there are two 
possible situations depending on N being even or odd. In the first case, there are 
two-minima that become deeper w.r.t other center minima. If N is odd, there is only 
one minimum. In general, the vacuum corresponds to non-trivial center vacua [7,12]. 
At low temperatures we recover the pure gauge potential. 

• Orbifold theories (c.f. [14]), for instance the Z2 case with SU(N) x SU(N) gauge 
group and fermions in the bifundamental representation (or theories with several 
gauge group factors in general) 

V(XW , A( 2 ) ) = £ V (A« - A«) + £ V (Af > - Xf >) - 2 J> (a« - xf >) . 

In the limit I — > 00 and large N, the diagonal sector A^ = A*- 2 -* becomes supersym- 
metric. At high temperatures the minima are on A^ = — A- 2 '' [12], and in this sector, 
the potential is equivalent to the orientifold case. 
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• Theories with more exotic matter content, etc. For a SU(N) gauge theory, whenever 
we have a representation of non-trivial iV-ality rjR, the center vacua will be shifted at 
high temperatures to a band with levels of r]R degeneracy and a splitting proportional 
to rm/N [7]. At low temperatures any contribution is exponentially supressed unless 
conformally coupled scalars are included. In that case, there will be extra zero modes 
that can modify the low-temperature regime. 



3. The G 2 case 

The results of section ^ can be applied to a gauge theory with exceptional group G2. We 
briefly summarize its relevant group theoretic properties (we refer to [15] for a more in-depth 
discussion of G2 gauge theory). G2 has rank 2, so the potential ( 2.16| ) is defined over a two- 



dimensional plane. The fundamental representation is 7 and the adjoint 14 dimensional. 
The group SU(3) is a subgroup of G2. The adjoint of G2 under this embedding splits into 
an adjoint and the fundamental and anti-fundamental representations of SU(3): 14 — > 8 
© 3 © 3*. Therefore, the roots of G2 are the same as the roots and the weights of the 
fundamental and anti-fundamental representation of SU(3). 

The fundamental representation of G2 can be decomposed in terms of SU(3) represen- 
tations as 7 — > 1 © 3 © 3*. Then, the Polyakov loop of G2 can be written in terms of a 
SU(3) matrix U 

( l \ 

L=tr U =l + 2RetrC7. (3.1) 

V u+ ) 

As explained in section ^, we choose to work in the gauge where the SU(3) matrix for flat 
connections is parameterized by the constant Cartan subalgebra of the group 

w x -6 Ju 2 -C Ju 3 -C 



U= diag [e w °, e lv % e lV 1 , (3.2) 

where the weights of the fundamental representation of SU(3), iT 1 , are not linearly inde- 
pendent but satisfy the following conditions: 



1 3 



3 

The physical region of the C plane is obtained quotienting by global periodic transfor- 
mations, that shift the phase of the eigenvalues of the Polyakov loop by multiples of 2ir: 
v 1 ■ C — > v 1 ■ C + 27rn*, and by Weyl transformations, that permute the eigenvalues of the 
Polyakov loop. In both cases, the root lattice determines the structure of the Cartan torus. 
Periodic transformations can be seen as translations along the roots a 13 = v l — v 3 , while 
Weyl transformations are seen as the reflections that move fundamental weights among 
themselves. The fixed lines of Weyl reflections lie along the roots. As a basis for the 
Cartan torus we can choose the simple roots a± = (l/y/2, ±^3/2). Since v % ■ a± = 0, 1, 
the period is 2ir. If we consider other representations, the group of Weyl reflections will be 
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different and the associated Weyl lines will show a different structure. For a SU (3) theory 
with fields in the adjoint representation we can also make Z3 center transformations that 
shift the eigenvalues by multiples of 2n/3. On the Cartan torus, this is seen as translations 
along the weights by 2ir units. In figures |l|, Q we show some examples. 




Figure 1: The Cartan torus (with roots normalized to unit norm) for SU(3) for the adjoint 
representation (on the left) and the fundamental representation (on the right). The Weyl lines are 
indicated by dashed lines. The global minima of the potential lie on the intersection of Weyl lines. 
For the fundamental representation, we also mark by white circles the smooth critical points of 
the potential. The center symmetry corresponds to translations on the torus that move the global 
minima of the adjoint representation among each other, notice that it is not a symmetry of the 
fundamental. 




Figure 2: The potential for G2 Yang-Mills at high temperatures. There is a global minima located 
at the origin of the Cartan torus and two local minima inside. 

The potential is then 

V{C) = 2Y,v(i?-c)+Y,v({t l -V)-c) . (3.4) 

In the SU(3) gauge theory, the first term in (|3,4j ) does not appear. In that case, the at- 
tractive contribution of the potential has its minima located where the cosine is maximum, 
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i.e. at 

(if - iP) ■ C = 2m , neZ. (3.5) 

For the adjoint representation, there are three minima located inside the physical region. 
In the pure Yang-Mills theory, one can move from one minimum to the other using Z3 
center transformations. Therefore, the localization around the minima at high tempera- 
tures could be interpreted as the weak coupling analog of the breaking of the center in the 
infinite volume theory. When the temperature is very low, the repulsive contribution of 
the potential dominates, so the minima become maxima and we expect the configurations 
to be spread along the Cartan torus, giving a vanishing expectation value to the Polyakov 
loop. When we consider matter in the fundamental representation of SU(3), or in the G2 
pure gauge theory, the center is no longer a symmetry. In the one- loop potential (3^), this 
is reflected by the presence of the first term. The minima of the attractive contribution 
are located at 

u 1 ■ C = 2im , n G Z . (3.6) 

So there is only one minimum inside the physical region. Then, only one of the minima 
of the SU{3) adjoint contribution is a global minimum, while the other two can only be 
local minima at most. However, the localization around the minimum still occurs when we 
change the temperature. 

In figure ^ we have plotted the G2 potential for different temperatures. We notice 
that in the high temperature phase there is a single global minimum located at the corners 
of the physical region. This minimum corresponds to the trivial vacuum, with trU = 1. 
There are also two local minima that in the pure SU(3) gauge theory will be the minima 
corresponding to the Polyakov loop values tr U = e ±27 ™/ 3 . When we lower the temperature, 
the minima become repulsive points due to the logarithmic term in the potential. The most 
repulsive point is the trivial vacuum, so the configurations will be distributed around the 
center values. 

We can compute the expectation value of the Polyakov loop at different temperatures. 
In terms of the eigenvalues Aj = Vi ■ C, the zero temperature T = path integral is 

Aj A t — r . 9 ( A 




2 

The expectation value of the Polyakov loop is exactly zero 



)ll sin2 (y)- 1 MX). (3.7) 



( l )t=o 



±-Jdfi(\) U + 2^cos(A i )j =0. (3i 



At non-zero temperature T = 1/(3, the Kaluza-Klein modes on S 3 give non-trivial contri- 
butions 

-21(1+2) 

. (3.9) 



wo 



1=1 



n 1+ 



sir ' 



, , v sinh* (^il) 




The expectation value becomes nonzero, and for high temperatures the measure is peaked 
around Aj = 0. This qualitative picture agrees with recent lattice results [16]. 
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(a) (b) (c) 

Figure 3: We show the Gi one-loop potential on the Cartan torus for (a) high temperatures 
(i/R = 0.01, (b) intermediate temperatures (3/R — 1.25 and (c)low temperatures (3/R = 100. 
As we lower the temperature, the initial single minimum Aq = becomes a maximum and new 
non-trivial symmetrically located minima appear, being pushed apart. The appearance of strong 
repulsion lines along the symmetry axis is due to the logarithmic term. These two features do not 
appear in other topologies, like the torus [7, 17]. 



4. Orthogonal and Symplectic Groups 



In this section, we will consider the phase structure of SO(2N) and Sp(2N) gauge theories. 
The center of Sp(2N) is Z 2 and that of SO(N) is Z 2 x Z 2 or Z 4 . for N even and Z 2 for N 
odd. However, since we are not considering spinors, only a %2 subgroup of the center will 
be relevant for our analysis for orthogonal groups for even TV. For odd TV, SO(N) groups 
have a trivial center. 

The non-trivial roots of SO(2N) and Sp(2N), in terms of the fundamental weights v, 
are given by ±z?j ± Vj (for SO(2N), i / j whereas there is no such restriction for Sp(2N)). 
For SO(2N), using the gauge transformation, the Polyakov loop can be brought in the 
form: 

/ cos#i sin#i \ 
— sin 6*i cos #i 

cos 6*2 sin 02 
— sin 2 cos 02 



O 



\ 



(4.1) 



COS0N sin0/v 
— sin#Ar cos0nJ 

Also, by a gauge transformation we can flip the signs of , s: 0{ — > —0{. The effective 
potential for the ^'s is given by 



V {0 l ) = Y,m-o 3 )+v{0 l +0 J ). 



(4.2) 



Introducing the density of eigenvalues p{0) = ^ + X^n^=i Pn cosn0, (where the Fourier 
coefficient p n is now real and p(0) = p(—0)), and converting summations to integrals, we 
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obtain 



V[p]=N 2 ^ d6p(6) ^ d9'p(9')f^-(l-z v (q n ))(cos(n(9-9')) + cos(n(9 + 9'))). (4.3) 

J-TT J-TT n=1 n 



Defining the Fourier components of p(9) as in (2.25), the potential can be written as 

' (l-z v (q n )) 



V[ Pn ] = iW £ ■ 



n=l 



n 



-pI- 



(4.4) 



For low temperatures, p n = 0, which implies a uniform distribution p(9) = and a 
vanishing expectation value of the Polyakov loop. At extremely high temperatures, the 
arguments of the cosine on the right hand side of (4.3) have to be minimized, which 
happens when the 9[ are the same and equal to or tt: i.e. 

p(0)=5(9) or 6(9 - tt). 

In the two saddle points, the Polyakov loop tr O = 2N J d9 p(9) cos(#) is 



trO = +1 
trO = -1 



p{9) = 5(9), 
p(0) = 8(9 -tt). 



If we sum over the different saddle points, we get a trivial value of the Polyakov loop. 
However, in the large N limit, we can choose one saddle point since large N provides a 
thermodynamic limit. In that case, the Polyakov loop is non-zero and the center symmetry 
is spontaneously broken. 

The situation is different for SO(2N + 1). The group SO(2N + 1) has no center. The 
roots, in terms of the fundamental weights z?, are given by ±i?j ± Vj and ±i?j. The Polyakov 
loop is now a (2iV + 1) x (2N + 1) matrix similar to ( [4.1[ ) with an extra row and column 
having 1 along the diagonal and zeros everywhere else. The effective potential for the #j's 
is now given by 

V(9i) = V(9i - 9j) + V{0i + 9j) + 2 V(9i) . (4.5) 

i^j i 

As before, we can introduce a density of eigenvalues and its Fourier components yielding 
an effective potential (4.4). The effective potential then is 

V[p] = N 2 f^\l - z v (q n ))( K J^ d9p(9) J W d9'p(9')(cos(n(9 - 9')) + cos(n(9 + 9'))) 



+ i f d9p(9) cos(n9) 

* " J — 7T 



(4.6) 



Defining the Fourier components of p(9) as in (2.25), the potential can be written as 



V\p n ) = N\ 2 Y, 



71=1 



:i-M<z")) (p2+ i 

n ttJM 



(4.7) 
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This is minimized when p n = —^n- This corresponds to an eigenvalue distribution p(9) = 
5^(1 + 27v) ~~ JN^i^)- The Polyakov loop is given by 

tiO = l + 2N f d6p(6)cos9 = l + 2N [ d9(^-(l + - ^5(9)) cos0 = 0. (4.8) 

J -TV J -TV 27r ^ 



In the extreme high temperature phase, because of the second line in ( |4.6| ) which picks 
out the 9 = saddle point, p(9) = 6(9), i.e. there is a unique saddle point in this case, 
giving a non-zero value of the Polyakov loop. The linear term on p n is similar to the one 
present for SU(N) theories with fundamental matter. In that case, the potential becomes 

V[p n \=Nir^ Pn+ ~N~n~ Pn ' ( ' 



n=l 



and is minimized at 



o ~ Nf ZF{qn) (4 10) 

Pn ~ 2Nl-z v (q n )' [ ' 

The value of the Fourier components is bounded to \p n \ < 1 because p(9) > 0. The 
bound is saturated at a temperature below the Hagedorn temperature, and the eigen- 
value density changes to a gapped distribution in a continous way. It can be shown that 
the deconfinement transition becomes Gross- Witten like [18,19]. However, in the case of 
SO(2N + 1) theories, both the quadratic and the linear coefficient come from the adjoint 
contribution and are the same. Hence p n = —l/(2irN) is fixed and therefore the transition 
is Hagedorn-like. 

In the large iV limit, the behavior of the theory with gauge group Sp(2N) is identical 
the the SO(2N) case. 



5. Gauss law and Screening 



In the previous two sections, we find that the Polyakov loop can be non-zero at high 
temperature is theories with no surviving center symmetry, even on finite volume. We are 
introducing a non-Abelian charge in finite volume and it does not cost infinite free energy. 
This is consistent with the non-Abelian Gauss law which we can write as 

j s (y i Et-f abc A\Ei)=Q a , (5.1) 

where Ef are the electric components of the non-Abelian field, is the magnetic potential, 
f abc are the group structure constants and Q a is the total external charge. The first term 
on the left hand side is a total derivative and vanish when integrated over a compact space. 
This is the only term that is present for Abelian fields and implies the known result that 
we cannot introduce an Abelian charge on finite volume. However, for non-Abelian fields 
the second term can be non-zero and Gauss law can be satisfied in the presence of a net 
external charge. 

For a gauge theory with group G and an unbroken center C, both the measure and the 
action are invariant under center transformations. This implies that there is a global center 
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symmetry which cannot be spontaneously broken in finite volume. Then, the expectation 
value of any center-breaking operator should be zero. The Polyakov loop in any representa- 
tion charged under the center transforms non-trivially under the center (which can be seen 
by performing large gauge transformations that are periodic up to an element of the cen- 
ter). Therefore, its expectation value should be zero, which corresponds to the statement 
that the free energy cost of introducing external particles that are charged under the center 
is infinite. When there is no center symmetry because the gauge group is centerless or the 
representation of the matter fields explicitly breaks it, there is no symmetry protecting the 
expectation value of the Polyakov loop and indeed we observe that it is nonzero in general. 

In the zero temperature limit, we find that the Polyakov loop is zero for any theory. 
We can also understand this in terms of Gauss law. In this limit, we can ignore the 
Polyakov loop and the term coming from the non-Abelian nature of the group in Gauss 
law. This is because we are essentially working in the zero coupling limit. The only non- 
trivial interaction terms we have kept correspond to the interaction of the fields with the 
background field a which is a zero mode of the theory. 4 There is no a at zero temperature 
and therefore no non-Abelian term in Gauss law. In this case, it costs infinite action to 
introduce any external charge and the Polyakov loop expectation value is zero. 

This is consistent with the physical picture in terms of screening. In the absence of a 
surviving center symmetry, an external particle in any representation can be screened by 
an appropriate number of gluon or matter particles, which allows us to introduce them at 
finite free energy cost. In theories with a center symmetry, an external particle charged 
under the center cannot be screened and therefore cannot have finite free energy in finite 
volume. 

6. Conclusions 

In this paper, we have investigated confinement on S 3 x S 1 for a gauge theory based on an 
arbitrary semi-simple Lie group containing matter in generic representations. Our analytic 
calculation reproduces the behavior of the Polyakov loop computed by lattice simulations 
for any theory for which lattice data exist. [16,17,20-22]. 

By using the general expression for the one-loop effective potential, we have shown that 
at zero temperature the Polyakov loop is exactly zero for any theory. Since the Polyakov 
loop is naturally related to the free energy of an isolated quark, our result implies that on 
a S 3 x S 1 manifold, any gauge theory confines. Although in this case confinement is due 
to kinematics, this finding is likely to be the weak coupling equivalent of confinement in 
gauge theories. For small but non-zero temperatures, the Polyakov loop is still zero if the 
theory has an unbroken center symmetry, but acquire a non-zero (albeit exponentially small 
0{e~P/ R )) vacuum expectation value otherwise. The very same features characterise the 
high temperature phase, although their explicit realisation differ from the low temperature 
regime. In the case in which the system is center symmetric, this symmetry cannot be 
broken on a finite volume, hence the Polyakov loop will still be zero at high temperature. 

4 In a normalization of fields that correspond to canonical kinetic terms in the action, this corresponds 
to taking the g — > limit after rescaling a — » -a. 
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However, this zero is due to a different phenomenon, namely the emergence of center- 
breaking minima of the effective potential, over which we have to take an average when 
we are not the in the thermodynamic limit. In this limit one minimum is selected, and 
this gives a non-zero value to the Polyakov loop. This is the case for SU(N) and SO(2N) 
gauge theories. For centerless theories, like G2 and SO(2N+ 1), there seems to be only one 
minimum for the effective potential in the high temperature phase. On this saddle point, 
the Polyakov loop is different from zero. 

In this paper, we have not considered Polyakov loop in the spinor representations 
of SO(N). Since these Polyakov loops will be non-trivially charged under the center of 
the gauge group, they should exhibit a behavior similar to the fundamental Polyakov 
loop in the SU(N) theory, i.e. at finite N, they should be zero at all temperatures. In 
the thermodynamic (large N) limit, after restriction to one center-breaking minimum, 
we should find a non-zero value of the Polyakov loop. We have also not studied the 
Polyakov loop in the adjoint representation of SU(N) theories. We expect to find that 
all representations uncharged under the center would have a non-zero Polyakov loop at 
non-zero temperature with a smooth cross over between the low temperature and the high 
temperature regimes. 

The minima of the effective potential exhibit a different structure in the low and the 
high temperature phase. This is true independently of the existence of an unbroken center 
group of the theory. This means that with a finite number of degrees of freedom, a cross 
over between the confined phase and the deconfined phase should take place. In the large 
N limit, this cross-over becomes a real phase transition. For finite N or in the case of 
exceptional groups, whether there is a cross-over or a real phase transition cannot be 
established by a computation on S 3 x S 1 . Recent lattice calculations suggest that in pure 
gauge theories the change of properties of the vacuum is always associated with a proper 
phase transition [17,20-22]. 

Another interesting question is whether calculations on a S 3 x S* 1 manifold can shed 
some light over the mechanism of color confinement. For QCD with matter in the adjoint 
representation, an interesting development in this direction is reported in [23]. 
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